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Abstract 

Although standard planar double bubbles are stable in the sense that 
the second variation of the perimeter functional is non-negative for all 
area-preserving perturbations the question arises whether they are dy¬ 
namically stable. By presenting connections between these two concepts 
of stability for double bubbles, we prove that standard planar double bub¬ 
bles are stable under the surface diffusion flow via the generalized principle 
of linearized stability in parabolic Holder spaces. 
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1 Introduction 

The standard double bubble is stable in the sense that the second variation 
of the area functional is non-negative. This follows for example from the fact 
that it is a local minimum of the area functional under volume constraints. It 
is however an open problem whether the standard double bubble is stable for 
volume conserving geometric flows such as the surface diffusion flow. 

The related problem for one bubble has been studied by Escher, Mayer 
and Simonett, see isizl, who showed that spheres are stable under the surface 
diffusion flow and the volume preserving mean curvature flow. In this paper we 
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show that the standard double bubble in is stable under the surface diffusion 
flow. In case of equal areas the result is illustrated in Figure [T] 

Before moving on to define the problem more precisely, let us make one 
point clear: Consider a (cost) functional having local minimizers. Even though 
minimizers exist it is not clear that an associated gradient flow will converge 
to these minimizers, see [5] for ODE examples. In other words, if a stationary 
state of the associated gradient flow is a local minimum, this in general does 
not imply stability of this equilibrium under the flow. 

As just mentioned, the surface diffusion flow is the volume preserving gradi¬ 
ent flow of the area functional. Indeed, it is the fastest way to decrease area while 
preserving the volume w.r.t. the id“^-inner product; see e.g. [Ulllin]. Let 
us now define the flow precisely. A surface is evolving in time under the surface 
diffusion flow if its normal velocity is equal to the negative surface Laplacian 
of its mean curvature at each point, that is, if a surface r(f) satisfies 


V{t) = -Ar(t)Lfr(t) • (l-l) 

Here V stands for the normal velocity, H is the mean curvature, and A is the 
Laplace-Beltrami operator, of the surface T(t). Surfaces with constant mean 
curvature are stationary solutions of the flow (HD. This flow leads to a fourth 
order parabolic partial differential equation (PDE). Thereby one may try to use 
PDE theories to answer the question on the stability of stationary solutions. 



Figure 1: An illustration of the stability of standard planar double bubbles, 
possibly up to isometries, (cf. the cover page to G. Prokert’s PhD thesis [15| 1 

Recently Priiss, Simonett and Zacher nsiini introduced a practical tool to 
show stability for evolution equations in infinite dimensional Banach spaces in 
cases where the linearization has a non-trivial kernel. It is called the generalized 
principle of linearized stability. This principle is extended in [1] to cover a more 
general setting. According to this principle, to prove stability, one needs to 
verify four assumptions known as the conditions of normal stability: 

(i) the set of stationary solutions creates locally a smooth manifold of finite 
dimension. 


2 









(ii) the tangent space of the manifold of stationary solutions is given by the 
null space of the linearized operator, 

(iii) the eigenvalue 0 of the linearized operator is semi-simple, 

(iv) apart from zero, the spectrum of the linearized operator lies in C+. 

We will see that the non-negativity of the second variation of the area functional 
plays an important role in verifying most of these assumptions for the double 
bubble problem. 

Let us note that the center manifold theory is used in [SI [7] to prove the 
stability of spheres under the surface diffusion flow and the volume preserving 
mean curvature flow. We remark that sofar no center manifold theory exists in 
the case of non-homogenous boundary conditions. Due to the triple junctions, 
we indeed get nonlinear boundary conditions in the corresponding PDE. 

Outline. In Section [5] we precisely dehne the problem which we summarize 
here: Let r° be an initial planar double bubble. We suppose that r° moves 
according to the surface diffusion flow including certain boundary conditions on 
the triple junctions. We continue then by observing that the set of stationary 
solutions consists precisely of all standard planar double bubbles. 

Next we transfer, via suitable parameterization, this geometric problem to 
a system of fully nonlinear and nonlocal partial differential equations with non¬ 
linear boundary conditions defined on fixed domains. We then linearize this 
nonlinear system. This is done in Section [3] 

In Section 15.11 we rewrite this nonlinear system as a perturbation of the 
linearized problem. We then see how suitably the problem fits to the generalized 
principle of linearized stability setting which is summarized in Section |4l 

It then remains to check the conditions of normal stability. Let us note 
here that understanding the geometric interpretations of the problem was of 
great help. Lemma lS. 131 proves assertion (iv). The non-negativity of the second 
variation is the main ingredient in the proof. Semi-simplicity is also proved by 
the non-negativity of the second variation in Section 1331 We prove assertion (i) 
in Section [5^31 and Corollarv l5.26l proves assertion (ii). 

By applying the generalized principle of linearized stability we then complete 
the proof of the stability, as summarized in Section |6l We continue then in Sec¬ 
tion IHT] to discuss general area preserving geometric flows. We then conjecture 
that the standard planar double bubbles are stable under sufficiently smooth 
area preserving gradient flows, see Conjecture 16.21 

In addition. Appendix El shows that the second variation is negative for two 
elements of the basis of the null space which correspond to non-area preserving 
perturbations. 
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2 The geometric setting 

A planar double bubble F C consists of three curves Fi, r 2 , Fa meeting two 
common points (triple junctions) at their boundaries such that Fi and 

F 2 (resp. F 2 and F 3 ) enclose the connected region Ri (resp. R 2 ). Hence the 
curve F 2 is the curve separating i?i and R 2 , see Figure [H 



Figure 2: A good example of a planar double bubble F = {Fi,F 2 ,F 3 } 


We study the following problem introduced by Garcke and Novick-Cohen 
m- Find evolving planar double bubbles F(t) = {Fi(f),F 2 (t),F 3 (t)} with the 
following properties: 


Vi = —Ar^Ki 

<(Fi(t),F 2 (t)) = <(F2(t),F3(t)) = <(F3(t),Fi(t)) 

Kl + K2 + K3 = 0 

VpiKi • noTi = Vr2«;2 • ndr2 = VrjKa • noTs 

F,(t)|t=o = r°, 


on Fi(t) 

^ on E(t), 
on S(t), > 
on S(t), 


( 2 . 1 ) 


where i = 1,2, 3, Fi(t) C and 

aFi(t) = dT2{t) = dTsit) ( = {p+it),p.{t)} =: m) • 

Here Vi is the normal velocity, Ki is the curvature, and Ar^ is the Laplace- 
Beltrami operator of the curve F^ {i = 1,2,3). Also Vr^ denotes the surface 
gradient and ngn denotes the outer unit conormal of F^ at dTi {i = 1, 2, 3). 

Moreover F° = {F]’,F 2 ,F 3 } is a given initial planar double bubble, which 
fulfills the angle ( 12 . 11 ) ^. the curvature (| 2 . 1 |l n and the balance of flux condition 
as above and satisfies the compatibility condition 

Apos:]^ Apo ^^2 F Apo Atg — 0 on H( 0 ). ( 2 . 2 ) 

Furthermore, the choice of unit normals ni{t) of Fi(t) is illustrated in Figure [31 
which in particular determines the sign of curvatures ki, K 2 and K 3 . We say that 
the curve has positive curvature if it is curved in the direction of the normal. 

Let us give a motivation for assuming the condition (12.21) on initial planar 
double bubble. 
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Figure 3: The choice of the normals 


Lemma 2.1. For a classical solution of the surface diffusion flow mi we have 

3 

^^AriKi = 0 onY,{f). (2-3) 

i=l 

Proof At the triple junctions p±(t) we can write for the normal velocities 




,ni{t)) . 

T—t 


Now the angle condition implies 




2=1 


^,n,(t)) = {— p ±{ t ) 


,^n*(t)) = 0. 


2=1 


As Vi = Ar^Ki, we obtain (lOD . □ 

Therefore if one seeks for a classical solution which is continuous up to the 
time t = 0, one should impose (12.311 on the initial data. 

After introducing the problem, let us see its interesting geometric properties: 

Lemma 2.2. A classical solution to the surface diffusion flow mi decreases 
the total length and preserves the enclosed areas. 

Proof. Assume r(t) is a solution to the flow (12.11) and let 



1 ds 


denote the total length. A transport theorem (see e.g. [3l Theorem 2.44]) gives: 
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=0 




-E / |Vr,(t)K*pds<0, ° (2.4) 

i=i 


where we used all the boundary conditions. Note that the sum of the normal 
boundary velocities I'dTi vanishes due to the angle condition, more precisely, 


3 

2=1 



3 

'^nori{t,p±{t)) = 0 . 
i=l 


Moreover, the integral over Y,{t) = {p+(t),p_(<)} should be understood as a 
sum over its elements. 

Next, let us prove that the enclosed areas are preserved: It is a standard 
fact that (see e.g. [121 equation (3.1)]) 


dt 


/ ldx= Vids — V 2 ds 

dri(t) Jr2(t) 


/ Ari(t)Kids+ / Ar2(t)K2ds 
/ Vr,(t)Ki • nor + / Vr,(t)K2 • ng = 0 . 

J'E(t) 


Similarly, we get ^ lR 2 {t) Idx = 0, which completes the proof. □ 

Let us mention that, via formally matched asymptotic expansions, the flow 
(|2.ip is derived as an singular limit of a system of degenerate Cahn-Hilliard 
equations in m, where in particular the boundary conditions at each triple 
junction are derived. 


2.1 Equilibria 

Let a planar double bubble L = {ri,r 2 ,r 3 } be a stationary solution of the flow 
(EH, i.e., r satisfies EH with Vi = 0 for i = 1, 2, 3. As a consequence 

Ar,K*=0 (i = l,2,3). 
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By the same arguments used in (12.41) we get 


0 = 


3 

E 

2=1 ' 


ds = - 


3 

E 

2=1 ' 


|Vi 


.KiPds. 


Thus VrjKz = 0 on r^. Therefore ki, /C 2 , K 3 are constant. Summing up, a planar 
double bubble T is a stationary solution of the flow (12.11) if and only if 

(i) the curvatures Ki are constant, with ki + K 2 K 3 = 0, and 

(ii) <(ri,rj) = ^ on E or equivalently X]i=i = 0 on S. 

It will turn out that the set of stationary solutions consists precisely of all 
standard planar double bubbles: 

Definition 2.3. A standard planar double bubble consists of three circular arcs 
meeting at their boundaries at 120 degree angles. (Here, we interpret a line 
segment as a circular arc too.) 

We refer to Figured for an example. Indeed, as circular arcs and line seg- 



Figure 4: The standard planar double bubble 

merits are the only curves with constant curvature, it just remains to verify the 
condition on curvatures. This is done in the following proposition given in [H 
Proposition 2.1]: 

Proposition 2.4. There is a unigue standard planar double bubble (up to rigid 
motions, i.e., translations and rotations) for given areas in The curvatures 
satisfy ki + K2 + K3 = 0 . 

Remark 2.5. As the choice of the normals in differs from ours, some sign 
differences particularly for the curvature quantities can occur. 

Therefore the set of all standard planar double bubbles 115 ^, 7 , 0 ( 01 , 02 ) forms 
a 5-parameter family (see Figure O, where 

(i) r > 0 is the radius of Fi, corresponding to scaling, 

(ii) ( 01 , 02 ) is the center of Fi, corresponding to translation. 
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Figure 5: The standard planar double bubble F = Z)i?r, 7 ,e(ai) 02 ) 


(hi) 


the angle 0 corresponds to counterclockwise rotation around the center of 

ri, 


(iv) the angle 0 < 7 < ^ corresponds to the curvature ratio. 

Indeed, by the law of sines we have for 77 ^^ 

Kl _ K2 _ K 3 
sin( 7 +^) sin (7 — |) sin (7 — tt) 


(2.5) 


and in case 7 = f we observe K 2 = 0 and ki = —K 3 . Note that due to our 
choice of normals we always have /ti < 0 and K 3 > 0. Moreover, 


K 2 > 0 for 7 < I, 
K 2 < 0 for 7 > I . 


For later use we define the constants qt as follows: 


q^ ■■= 



Kfc) 


for (z, j, k) = (1, 2, 3), (2, 3,1) and (3,1, 2). Then the following result is true. 
Lemma 2.6. We have 


rcot(7-f)/C2 7 7^f, 

= C0t(7 + |)ki, q 2 = l ^ 93 = cot(7 - 7r)K3 . 

iiRI) ^ = 
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Proof. We calculate 



^ 2 / sin( 7 +|)cos(§) 
73 V sin( 7 -f) 



and obviously q 2 = ^ki = for 7 = |. 

formula (|2.5I) . The proof for <71 and <73 is similar. 


The continuity follows from 


□ 


Moreover, using the sum-to-product trigonometric identity, we get 



( 2 . 6 ) 


One strategy to deal with geometric flows on hypersurfaces is to parameterize 
the evolving hypersurfaces with respect to a fixed reference hypersurface. This 
eventually leads to a PDE on a fixed domain allowing us to employ PDE theories. 

3 PDE formulation and linearization 

In this section we introduce the proper setting to reformulate the geometric 
flow ( 12 . 1 |) as a system of partial differential equations for unknown functions 
defined on fixed domains. For this, we employ a parameterization with two 
parameters. The parameters correspond to a movement in normal and tangen¬ 
tial directions. This parameterization is adapted for two triple junctions from 
Depner and Garcke [4], see also [5]. 

3.1 Parameterization of planar double bubbles 

Let us describe ri(t) as a graph over some fixed stationary solution P* using 
functions 


p, :r* X [0,r)^R (* = 1,2,3). 

The precise way how pi defines Ti(t) will be derived in what follows. 
Fix any stationary solution 


P* = DBr*^j^^ 04 al,a* 2 ) (r* > 0 , (a^a^) 0 < 7 * < ^, 0 < 0 * < 2 **) . 


Then we observe 
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where 21* is the length of F* {i = 1,2,3) and of course k* = — 

Let <!)* : [—l*,l*] —>■ be an arc-length parameterization of F*. Hence 

F* = {$*(a;) :x€ 

Furthermore, set = x{a) G R, for a G F*. To simplify the presenta¬ 

tion, we hereafter set 

daw{a) := dx{wo^*){x), a = ^*{x), (3.1) 

that is, we do not state the parameterization explicitly. Also we slightly abuse 
notation and write 

w{a)=w{x) (ctGF*). (3.2) 

To parameterize a curve nearby F*, define 

T, : F* X (-e, e) X (-5, 5) ^ R2 , (3.3) 

(cr, w, r) I—'I'i(CT, w,r) := a + w n* (cr) + r r* (u). 

Here t* denotes a tangential vector field on F* having support in a neighborhood 
of 9F*, which is equal to the outer unit conormal nor* at 9F*. 

Define then = ($i)pt w often omit for shortness the subscript (pi, /x^)) 

by 


: F* X [0,r)R^ ■.= ^,{a,p,{a,t),^i,{w^{o),t)), (3.4) 

for the functions 

p, :F* X [0,T)^ (-e,e), p, : E* x [0, T) ^ (-5, <5), (3.5) 

where, similarly as before, S* = 9F* = {p*j^,p*_}- 

The projection pr^ : F* —>■ E* is defined by imposing the following condition: 
The point prj(cr) G 9F* has the shortest distance on F* to a. Clearly, in a small 
neighborhood of 9F*, the projection pr^ is well-defined and this is sufficient for 
us since this projection is just used in the product t)T* (a), where the 

second term vanishes outside a (small) neighborhood of 9F*. 

Now let us set, for small e,d > 0 and fixed t, 

($,)t : F* ^ R2, ($0t('^) := t) Vu G F* 

to finally define a new curve 

rpi.w(i) := image(($*)t). (3.6) 

Observe that for pi = 0 and pi = 0, the curve Fp._^. (t) coincides with F* for all 

t. 

At each triple junction, we have prepared for a movement in normal and 
tangential direction, allowing for an evolution of the triple junctions. Therefore, 
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we can now formulate the condition, that the curves ri(<) meet at the triple 
junctions at their boundary by 


= ^ 3 {cr,t) for cr e S*, t > 0 . (3.7) 

Next we prove that this condition leads to a linear dependency at the boundary 
points. As a result, nonlocal terms will eventually enter into PDE-formulations 
of the geometric evolution problem. 

Lemma 3.1. Equivalent to the equations (1X71) are the following conditions 
( (i) 0 = pi+P 2 + P 3 on S*, 

I t"\ ^ t \ V* (^■^) 

( (ii) h-i =-^{Pj - Pk) on S , 

for {i,j,k) = (1,2, 3), (2, 3,1) and (3,1,2). 

Here the linear dependency (ii) can be recast as the matrix equation 

fi = Jp onE*, (3.9) 

with the notations p = (p-i, ^ 2 ,^ 3 )^ P = {pi, P 2 , P 3 ) and the matrix 


J 



1 

0 

-1 



Proof First we prove that (13.71) implies (13.81) . Using the definition of <i>i, (13.7p 
can be rewritten as 


Pi n* + Pi nor* = Pj n* + pj nor* on E* (3.10) 

for (i,j) = (1,2), (2,3). By setting 

q := Pi n*i + pi nar* =^ 2 ^ 3 + P 2 ngr* = P3n*3+ ps nor* on E* 
we obtain pi = (g, n*) for i = 1,2, 3. Thus the angle condition for T* gives 

'^P^ = ® • 

i—1 i—1 i—1 

This proves (i). As a result of (13.101) we see further 

Pi{n*,n*) + pifagri^n*) = pj on E* . 

On the other hand the angle condition implies 

{n*,n*) = cos{^), (nar*,'oP = cos(27r-(^ + |)) =-sin(^) on E* 

for {i,j) = (1, 2), (2,3), (3,1). Therefore using (i) we conclude 

Pi = -^Pj - cpf) = -^((1 + c)pj + cpk) = f (pj - Pk) , 

where s := sin(^) and c := cos(^) = —^ and this yields assertion (ii). The 
proof of the converse statement is explicitly given in [U Lemma 2.3]. □ 
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Note that we followed m in proving statement (i), while an easier proof is 
given here for assertion (ii). Notice further that (13.81) easily implies 

Ml + + M3 = 0 on E* . (3.11) 

Remark 3.2. Let us now note that it is within this set, i.e., the set of all planar 
double bubbles which can be described as the graph over F*, that we will seek a 
solution to the problem 

Naturally, we assume also that the initial double bubble r° from (12.ip is 
given as a graph over F*, i.e., 

r° = {«'*(cr,/9°(cr),M°(pr(CT))) : cr £ F*} (i = 1,2,3) 

for some function Here fjP = LT on E* as F° is assumed to be a double 
bubble, i.e., the curves F° meet two triple junctions at their boundaries. 

3.2 Nonlocal, nonlinear parabolic boundary-value PDE 

The idea is to first derive evolution equations for pi and Mi which have to hold 
if the Fi (z = 1,2,3) in (13.61) satisfy the condition (13.71) and solve the surface 
diffusion flow (EH) and then to make use of the linear dependency (|3. H p m 
deriving evolution equations solely for the functions pi. 

As you may have noticed, nonlocal terms will appear in the formulations 
since this linear dependency p3.9l) just holds at the boundary points. 

Appendix provides for the reader’s convenience the derivation in detail. 
Indeed a similar derivation is done in [T], which is originally given in 0- 
Therefore, let us present the final system of fourth-order nonlinear, nonlocal 
PDEs for t > 0, z = 1, 2, 3 and j = 1, 2,..., 6: 

{ dtPi = ^i{pi,p\T.-) 

+ ^^{pi,p\T.-){{J{I-‘^{ p,p\t.-)J) ^5’(m,p|e*)} oPH)i onF*, 

0 = 0j (m) on E*, 

(3.12) 

with the initial conditions 


Mi(-,0) = p° oiiF*, 

where in particular 'SiiPi, p\y.*) is a fourth-order nonlinear equation in pi. 

Remark 3.3. Note that the price to pay for obtaining eguations solely for func¬ 
tions Pi is the appearance of nonlocal terms, in particular the nonlocal terms of 
highest-order (fourth-order) 1 ?(m, pis*) o pri into the formulation. 

As demonstrated at the beginning of Appendix lEl the functions 3’i,f8i,0j 
are rational functions in the p-dependent variables, with nonzero denominators 
in some neighborhood of p = 0 in (F*) (can be inside of square roots equalling 
to 1 in some neighborhood of p = 0 in (^^(F*), see the term d-). 
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3.3 Linearization around the stationary solution 

The linearization of the surface diffusion equations and the angle conditions 
around the stationary solution p = 0 are done in [H Lemma 3.2] and [U Lemma 
3.4] respectively. 

Remark 3.4. Note that the situation in S' is slightly different from ours, but 
nevertheless the results obtained there are applicable to our problem. More pre¬ 
cisely, the authors in consider the situation where, apart from the appearance 
of a triple junction, one has to deal with a fixed boundary. However, as they 
assume that the triple junction will not touch the outer fixed boundary, they can 
use an explicit parameterization, exactly as ours, around a triple junction and 
another parameterization near the fixed boundary and finally they compose them 
with the help of a cut-off function. Thus we can use their result for each triple 
junction. 

Therefore, taking into account the linear dependency (ii) from Lemma 13.11 
we get for the linearization of the nonlinear problem (13.1211 around p = 0 (that 
is, around the stationary solution T*) the following linear system for i = 1, 2, 3 

^tp^ + Ar* (Ar.p, + {K*fp^) =0 in T* , (3.13) 

with the boundary conditions on E* 

Pi + P2 + P3 = 0 , 

Qi Pi + Pi ~ Pi "b Pj (b j) = (1; 2), (2, 3), 

= 0 , 

+ (^i) Pi) = ^rtgp. (Ar*Pj + (nj) Pj) ihj) = (1)2), (2,3), 

(3.14) 

where 

for (*,j,fc) = (l,2,3),(2,3,l) and (3,1,2). 

Let us recall the parameterization (remember our abuse of notation (13.21) 1 
and employ the following facts 

^r* Pi = dlp^ for X G [-1* ,1*], 

dngp,Pi = Vr*Pi • nor* 

= dxPi (T* ■ no ^,) = ±da:Pi at X = ±1* , 
l^rigr. = K* at X = ±1* , 

where x is the arc length parameter of T* and denote by T* the tangential 
vector of T*. We can then rewrite the linearized problem in terms of functions 
Pi ■ \~^i ) ^ [1^) 

dtPi + dl{dl + {K*f)p^ =0 for a; G [-l*,l*] 


13 






with the boundary conditions 

Pi + /02 + P3 = 0 , 

qlpi ± dxPi = q2P2 ± dxP2 = q^Pd. ± d^ps , 

^ Eli(5^P. + K*)V0=O, 

dx{dl + {K\f)pi = d:,{dl + {K*2f)p2 = d:,{dl + {nlf)p:i ■ 

In the boundary conditions (13.151) we have omitted the terms ±Z* in the functions 
Pi. That is, for instance the boundary condition pi +P 2 +P 3 = 0 should be read 
as 

Pi(±^i) + P2{±ll) + P3i±ll) = 0 . 

Furthermore, notice that the linearized problem is completely local as, in par¬ 
ticular, we linearized around a stationary solution. 

Now we are in a position to look for a suitable PDF theory in order to 
answer the question of stability. The generalized principle of linearized stability 
in parabolic Holder spaces, proved in |T], see also [Ml HU, provides the tool. 

4 The generalized principle of linearized stability 
in parabolic Holder spaces 

In this section we present the practical tool, proved in [T], for proving the sta¬ 
bility of equilibria of fully nonlinear parabolic systems with nonlinear bound¬ 
ary conditions in situations where the set of stationary solutions creates a C^- 
manifold of finite dimension which is normally stable. The parabolic Holder 
spaces are used as function spaces. 

Let fl C R" be a bounded domain of class ( 72 m-i-a mSN, 0<a<l 
with the boundary dVl . Consider the nonlinear boundary value problem 

{ dtu{t,x)A{u{t,-)){x) = F{u{t, .)){x), X G it, t>0, 

Bj{u{t,-))[x) = Gj{u{t,.)){x), xGdn, j = l,...,mN, ( 4 . 1 ) 

u(0, x) = uo(x), X € fl, 

where u : H x [0, 00 ) —>■ R^. Here A denotes a linear 2mth-order partial differ¬ 
ential operator having the form 

(Au)(x) = a.y{x)'V~*u{x), xGfl, 

|7|<2m 

and Bj denote linear partial differential operators of order rrij, i.e., 

{Bju){x) = ^ bjj(x)V^u(x), X G dfl, j = l,...,mA^, 
l/3|<mj 

with 0 < TOi < 1712 < • • • < 'm-mN < 2m — 1. 

The coefficients a-y{x) € R^^^, F^(x') € R^. We assume that 
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(H2) the elements of the matrix aj(x) belong to C°‘(fl) and 

the elements of the matrix 6^(x) belong to (917). 

Concerning the fully nonlinear terms F and Gj, let us suppose 

(HI) F : B{0,R) C C'^"*(17) —>■ (7(17) is with Lipschitz continuous deriva¬ 
tive, F(0) = 0, F'{0) = 0, and the restriction of F to B(0, R) C (7^'"“''“(17) 
has values in (7“(17) and is continuously differentiable, 

Gj : B{0, R) C (7™-» (17) —?> C(dfl) is with Lipschitz continuous second- 
order derivative, Gj(0) = 0, (7'(0) = 0, and the restriction of Gj to 
B{0,R) C (72’"+“(17) has values in (72™+“-'"^ (917) and is continuously 
differentiable. 

We set B = (Hi,..., H^at) and G = (Gi,..., G^w)- 

We denote by ,7 C Bxi{0, R) the set of stationary solutions of (j4.1L i.e., 

uG£ <;=^ u G Bxi{0, R), Au = F{u) in 17, Bu = G{u) on 917 , (4.2) 

where Xi = G^™+“(17). The assumption (HI) in particular implies that 

M* = 0 belongs to £. 

Now the key assumption is that near u* = 0 the set of equilibria £ creates a finite 
dimensional G^-manifold. In other words we assume: There is a neighborhood 
C/ C of 0 G (7, and a G^-function 4' : U Xi, such that 

• 4'([/) C £ and 4'(0) = it* = 0, 

• the rank of 'I''(0) equals k. 


Moreover, we at last require that there are no other stationary solutions near 
u* = 0 in Xi than those given by 'I'((7). That is we assume for some n > 0, 

£ABxAu*,ri) = ^{U). 

The linearization of (HU) at u* = 0 is given by the operator Aq which is the 
realization of A with homogeneous boundary conditions in X = G(17), i.e., the 
operator with domain 

I7(Ho) = juG G(17)n n IH^'"’^’(17) : Au G X, Hu = 0 on 917|, 

l<P<+oo 

Aqu = Au, u G D{Aq) . 


Let v{x) denote the outer normal of 917 at x G 917. We assume further the 
normality condition: 


for each x G 917, the matrix j 

{x){v{x)Y J 

^ where { i = 1 ,..., nj, } = { j : 775 , 5 = k } , 


is surjective. 


(4.4) 






Suppose at last the following first-order compatibility conditions holds: For 
j such that rrij = 0 and x S dVl 

I Bu° = Giu°), 

\ Bj{Au° - F{u°)) = G'j{u°){Au° - F{u°)). 

Theorem 4.1. ([1] Theorem 3.1]) Let u* = 0 be a stationary solution of (]4.1I) . 
Assume that the regularity conditions (HI), (H2), the Lopatinskii-Shapiro con¬ 
dition, the strong paraholicity and finally the normality condition ((Oil hold. 
Moreover assume that u* is normally stable, i.e., suppose that 

(i) near u* the set of eguilibria £ is a G^ -manifold in Xi of dimension k € N, 

(ii) the tangent space of £ at u* is given by N (Aq), 

(hi) the eigenvalue 0 of Aq is semi-simple, i.e., i?(Ao) ® N (Aq) = X, 

(iv) (T (^o) \ {0} C C+ = {z € <C : Re z > 0}. 

Then the stationary solution u* is stable in Xi. Moreover, if is sufficiently 
close to u* in Xi and satisfies the compatibility conditions (lOl) . then the prob¬ 
lem (HU) has a unique solution in the parabolic Holder spaces, i.e., 

we (7^+^. 2 ™+“([0,00) xH) 

and approaches some G £ exponentially fast in Xi as t —>■ 00 . 

Remark 4.2. We refer to JT] Section 2] for the definitions of Lopatinskii- 
Shapiro condition and the strong parabolicity as well as for a complete treat¬ 
ment. 

In order to apply this theorem to prove stability, we must first show that our 
nonlinear, nonlocal problem (13.121) has the form (14.11) . We then devote the rest 
of the paper to show that the problem p.l2|) verifies all hypothesis of Theorem 

mu 


5 Verifying the hypotheses of Theorem 14.11 


5.1 General setting 

If we change the variables by setting for each i = 2,3 


X = — I* 


21 * 


21 * 


xG [-11,11], 


then we easily can restate the nonlinear, nonlocal system (13.121) as a perturbation 
of a linearized problem, that is of the form mu, with £t = [—1*, I*], 


dip, 



■(li)" 

0 0 


0 0 

Ap = 

0 

ihY 0 

dtp + 

0 0 


0 

0 {h)\ 


0 0 ( 13 k: 5 ) 2 _ 
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and 


Bip = [1 1 

B 2 P = ± [li 

1 ] p, 

-I2 Q]d^p+ [g* -52 

0 ]p, 


H3P = ±[0 

I2 -I3] i9£>:P +[0 52 

- 93 ] p> 


BaP= [(li)" 

( 12 )^ {hf] dip + [{nlf 

(^ 2 )^ («3)^]p. 


B,p= [(li)3 

-(12)3 0 ] dip+[l,i 

F!)" -l2('^2)^ 

0 ]d^p, 

Bqp = [ 0 

(12)3 -(l3)3]a3p+[ 

0 12 (^ 2 )^ 

-l3(«^3)^] d^p 


To simplify the presentation, we have dropped the tilde. Here 


P ■ X [0,oo) 


P = [pi,P2,Pzf 


P 

and the constants are given as 1^ := 77 (* = 1, 2, 3). 

i 

When we write (j3.12|) in the form of (14.11) . the corresponding F is a smooth 
function defined in some neighborhood of 0 in C^[ri) having values in (7(11). The 
reason is that, F is Frechet-differentiable of arbitrary order in some neighbor¬ 
hood of 0 (using the differentiability of composition operators, see e.g. Theorem 
1 and 2 of m Section 5.5.3]). The same argument works for the corresponding 
functions Gi. We have obtained that assumption (HI) is satisfied. 

Obviously, the operators A and Bj satisfy the smoothness assumption (H2) 
and the operator A is strongly parabolic. Now Let us check that the Lopatinskii- 
Shapiro condition (LS) holds. To verify this, for A G C+, A 7 ^ 0, we consider the 
following ODE 


' Xv,{y) + {\ifd^v,{y) =Q, (y > 0 ), 
ui( 0 ) -I- ^2(0) -I- ^3(0) = 0 , 

li( 9 yz;i( 0 ) = \2dyV2{0) = l3i9j;U3(0), (5.1) 

= 0 , 

(ll)^(9^Ul(0) = (l2)^9yt'2(0) = (l3)^l9yU3(0) 


and we show that u = 0 is the only classical solution that vanishes at infinity. 
The energy methods provide a simple proof: We test the first line of the equation 
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(|5.1|l with the function —Vi and sum for i = 1, 2, 3 to find 



3 «oo 

^(h )3 / ^H^yV^dy 
„-_i Jo 


2=1 

3 


=0 


/ dyVidyVidy + '^v^ (li)^dyVi 


2=1 


2=1 


=0 


oo 

0 


w PUCJ ^ 

2—1 *^0 2 — 1 

3 «oo 

do 


\idyVi 


oo 

0 


Here we have used all boundary conditions at j/ = 0 and the fact that the func¬ 
tions Vi and consequently all their derivatives vanish exponentially at infinity. 
The latter holds due to the fact that the solutions of the above equations are 
linear combinations of exponential functions. The facts that 0 ^ A G C+ and 
li > 0 enforce u = 0. This verifies the claim. 

Furthermore, the matrices 


[1 1 1]> 


li —h 0 

0 I2 -I3 ’ 


[(li)" (I2)" (13)^] , 


(ii)^ -(b)^ 0 

, 0 (b)^ -(13)3 


are surjective and hence the normality condition (14.41) is satisfied. 


5.1.1 Compatibility condition 

We next turn our attention to the corresponding compatibility condition (14.51) . 
As we have assumed the initial planar double bubble T^ fulfills the contact, 
angle, the curvature and the balance of flux condition, we see yP = Pf and 
&j{p^) =0 for j = 1, 2,..., 6. This is exactly the first condition in (14.51) . 

Concerning the second equation in the compatibility condition (14.51) . the 
following lemma shows that it is equivalent to the geometric compatibility con¬ 
dition (j2.2p if the existence of triple junctions and the angle condition for the 
initial data are already assumed. 

Lemma 5.1. Under the conditions <3j{p^) = 0 {j = 1,2,3) and pP = Jp^ 
on E*, the second equation in the corresponding compatibility condition (H3]) 
and the geometric compatibility condition p.2l) are equivalent, provided p^ is 
sufficiently small in the C^-norm. 
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Proof. The second equation in the corresponding first-order compatibility con¬ 
dition (|4.5p reads as 

z: — 

E 'il-B{p°,p°)j)-^:S{p°,P°) ), = 0 (5.2) 

i=l 


on S*. Here we have used the facts that the zeroth-order boundary operator 
Biu = X]i=i G*! = 0. Let us remind that 








(<>»^°) 


and r* = on S*. 


On the other hand, the angle condition implies 

{nhni) = {n*,n°) , (ngp.,n°) = (n^p.,n°) on S*. 
Thus (15.21) can be rewritten as 


—^ Y,A{.,p<^, {Jp%)n. (., p°, {Jp\) + <Bi Y.{Jz)i = 0 on E* , 

i=i 

where {n\,n\) 7 ^ 0 if r° is close enough to T* in C^-norm, that is if is 
sufficiently small in the C^-norm. 

Moreover, due to the definition of the matrix we have 


3 

= o V2/eR3_ 

Hence the compatibility condition (|5.2[) is equivalent to 


i=l 

which is exactly the geometric compatibility condition (12.21) written in a param¬ 
eterization. This finishes the proof. □ 


5.2 The spectrum of the linearized problem 

Since fl = [—C M, the linearized operator Aq (see 63) is defined as 
Aqu = Au with domain 

D{Ao) = G C'‘*(H) : Bu = 0 on 9 h|, 

where A and B is defined in Section [5Tl Due to Remark 2.2 in [1], the spectrum 
of the linearized operator Aq consists entirely of eigenvalues. As the analysis of 
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the eigenvalue problem is invariant under the change of variables, we switch to 
the setting where the functions Ui (i = 1, 2, 3) have different domains. 

Now, the eigenvalue problem for the linearized operator Aq reads as follows: 
For i = 1, 2, 3, 

Ar* (Ar*Ui + {K*Yui) = Xui in F* {i = 1, 2, 3), (5.3) 

subject to the boundary conditions on S* 

Ui + U2 + U3 = 0 , 

Qi ^3 ^^ar’^ ^3 

< s o ' (5-4) 

+ {K*Yui = 0 , 

a„ar. (Ar*wi + {K*fui) = (Ar-Uj + {K*fuj) , 

V t 3 

where (i,j) = (1,2), (2,3). 

To derive a bilinear form associated with this eigenvalue problem, let us 
multiply the equation (|5.3p by — (Ar*ul+ {K*)'^ut) and then integrate by parts 
and sum over i = 1, 2, 3 to find 


3 » ^ r 

/ |Vr*(Ar*Ui + (K*)^Mj)|^ds = / Mj(Ar*u7+(O^wl) ds . 

i=i dr- 

Here, as usual, we have used the last two boundary conditions. We observe 
further 

/ ■Ui(Ar*u7 +(O^wT) ds = / |Vr-Mi|^ - (O^lu^p ds 

i=i dr- dr- 

^ ^ r 

^ ^ / '^i djig-p* UiUi . 

ds- * 


On the other hand 

3 


i=l 


'E- 


'y ' / Uidngj^tUi — ^ ^ / (ui l9ngp» Mi + IWil 9i|Mi| ) 
ds- • ds- 

3 „ 3 p 

= ^ ' / (^ngp» Ui + Qilli^Ui — 'y ( 

i=i ds* * 

P 3 3 

(^ngp. Ml + g>l) y] Ml - y] 

= -y^/ 








2=1 2=1 
=0 


2=1 


/s* 
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We now combine the three equalities above to discover 



3 


(5.5) 


where 



3 


3 


(5.6) 

Note carefully that in (15.61) we just used integration by parts to obtain the 


second equality. It is interesting now to see that although (due to the linearized 
angle condition and the fact that on the boundary + M 2 + W 3 = 0 ) we have 



(5.7) 


but nevertheless this does not effect the value of I{u,u) (cf. [T^ Remark 3.7]). 

Remark 5.2. The identity (15.51) in particular shows that A S R. 

Remark 5.3. Indeed as one may have expected, the linearized problem (|3.13|) . 
(|3.14|1 is the gradient flow of the energy functional 



with respect to the H ^-inner product, see for instance m- 

5 . 2.1 Related problem: Double bubble conjecture 

The goal of this section is to prove that, a part from zero, the spectrum of the 
linearized problem lies in ]R_|_. We do this by considering the bilinear form /(, ). 

In the following we state the second variation formula proved in general 
dimension by Morgan and co-authors: 

Proposition 5.4. 1 |121 Proposition 3.3]). Let P* be a stationary planar double 
bubble and let ipt be a one-parameter variation which preserves the areas of 
enclosed regions. Furthermore denote by L{t) the length of (ptfL*)- Then 



where m = 
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Here and hereafter, by (one-parameter) variations {^Pt}\t\<<L : T —^ of a 
double bubble F C K.^ we mean the variations which are smooth (up to the 
triple junctions) having equal values along triple junctions. 

Remark 5.5. Notice that in (EH) we have used outer unit conormals where 
inner unit conormals are used in m- In addition, the constants q* and their 
corresponding ones in m are also opposite in signs due to the different choice 
of normals. This explains the sign differences. 

Remark 5.6. Of course, a double bubble is stationary for any variation pre¬ 
serving the area of the enclosed regions if and only if it is stationary for the 
surface diffusion flow see Section \2.1\ and page f65]. 

Following [12], we denote by J^(r) the space of functions u G H^{T) satisfying 


ui -f M 2 + 1*3 = 0 on E , 


Ui= U2= U3 . 


/ri 






Lemma 5.7. 1 |12L Lemma 3.2]). Let F* be a stationary double bubble. Then 
for any smooth u G J^(F*) there is an area preserving variation {ift} o/F* such 
that the normal components of the associated infinitesimal vector field are the 
functions Ui, i.e., Ui = {-^(ptyn*), i = 1,2, 3. 

We are now ready to present: 


Definition 5.8 (The concept of stability in differential geometry). A double 
bubble F* is said to be variationally stable if it is stationary and 


l{u,u)>0 VuGJ^(F*). 


We are forced here to name the concept of stability in differential geometry 
variationally stable instead of stable. Indeed it is an open problem whether for 
double bubbles this concept of stability in differential geometry is equivalent 
to the concept of stability in PDF theory. There are several evidences in this 
work which show how closely these two concepts are, starting from Lemma [5.131 
below. 


Remark 5.9. Note that the concept of stability in differential geometry was 
called stable in Hi- 

Corollary 5.10. A perimeter-minimizing double bubble for prescribed areas is 
variationally stable. 

Proof. Let F be a primeter-minimizing double bubble. As a minimizer, the 
second derivative of length is nonnegative along all variations which preserve 
the area, in other words by Proposition 15.41 liu. u) > 0 for all functions u 
given by normal components of area preserving variations. On the other hand, 
by Lemma 15.71 we know that every smooth element of J^(F) is of this form. 
Therefore I{u,u) > 0 for all u G J^(F), which finishes the proof. □ 
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Theorem 5.11. ([51 Theorem 2.9]). The standard planar double bubble is the 
unique perimeter-minimizing double bubble enclosing and separating two given 
regions of prescribed areas. 

Therefore as an important corollary one gets: (see also [141 Theorem 3.2]) 
Corollary 5.12. The standard planar double bubble is variationally stable. 

We are now ready to see the first evidence. 

Lemma 5.13. (t(^o) \ {0} C M+. 

Proof. Let A S ^(^lo) \ {0}- As mentioned before the spectrum consists entirely 
of eigenvalues. In addition, according to Remark 15.21 A is real. 

Therefore, let A be an eigenvalue with a corresponding eigenvector u S 
( 74 +“(r*). This means u solves the eigenvalue problem (15. 3|) subject to the 
boundary conditions (15.41) for A. Since A 7 ^ 0, we deduce after integrating (15.31) : 



where we employed the divergence theorem and the last boundary condition. 
This together with the first boundary condition implies that u G J^(r*). There¬ 
fore I{u,u) > 0 by Corollarv l5.12l 

Now assume I(u,u) = 0. In view of the equation (15.51) . we obtain 

Ar*Mj -I- {K*iYui = Ci 

for some constants Ci (cf. [T2l Lemma 3.8]). This together with the equa¬ 
tion (|5.3|) immediately implies u G N{Aq), i.e., A = 0, a contradiction. Thus 
I{u,u) > 0 for the eigenvector u. Now A > 0 by (15.3|) . This finishes the 
proof. □ 

The bilinear form /(, ) is further discussed in Appendix [A] 

5.3 Null space of the linearized problem 

We next determine the null space of the linearized operator Aq. That is, we 
consider the case A = 0 in the eigenvalue problem (I5.3I) . (I5.4I) . 

Using the identity (15.51) . we easily get u G N{Ao) if and only if there exists 
a constant vector c = ( 01 , 02 , 03 ) G such that 

Ar*Ui{K*)^Ui = a on r* (i = 1,2,3), (5.8) 

subject to the conditions 

{ ui -I- U 2 + U 3 = 0 on E*, 

q^Ui -f dng^,Ui = q 2 U 2 -\- dng^,U 2 = q^u^i + dng^,U3 On E*, ( 5 . 9 ) 

Cl -I- 02 -I- 03 = 0 . 

Notice that the constant vector c = c(u) depends linearly on u by (15.81) . 
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Definition 5.14. Following JT^, we define the space of Jacobi functions 

J{T*) := {u G N{Aq) : c = c{u) = 0 } . 

We need, for later use, an identity that relates the null space N{Aq) to the 
bilinear form /(, ). 

Lemma 5.15. Assume u G N{Aq). Then 

^ ^ r 

I{U,U) = - / Ui, 

i=i 

where the constants Ci, satisfying = 0, depend linearly on u by (|5.8I) . 

Proof. By inserting (15.8p into the definition of the bilinear form (15.611 and taking 
into account the equation (15.71) coming from the first two boundary conditions 
in (15.91) . we get the desired identity. □ 

As a corollary we get 

Corollary 5.16. If u G N{Aq) (d then I{u,u) = 0. 

Let us rewrite the linear equations (15.81) as a system of linear nonhomoge- 
neous second order ordinary differential equations with constant coefficients 

dlui +{n'f'Ui = a ioT X e [-l*,l*] (* = 1 , 2 , 3 ), 

with the conditions 

{ Ml + 7*2 + W3 = 0 on E*, 

g*Mi ± dooUi = qlu2 ± d^U2 = q^us ± on E*, 

Cl + C2 + C3 = 0 

for the functions m : [—11, ?*] —>■ M. 

5.3.1 Determination of Jacobi functions 

Let us first consider the case 7^ 0- The general solution of the linearized 
problem is then 


Ui{x) = aiSm{K*x) + biCOs{K*x) (* = 1 , 2 , 3 ). ( 5 . 10 ) 


We calculate at a; = ±Z* 

qlui = =F cot(7* + |)Kiai sin(7* + f) + cot(7* + ^)Klbi cos(7* + |) 
= Toi^i cos(7* + f) + biK^ cot(7* + |) cos(7* + |), 

±dxUi = ±aiK,\ cos(7* + f) + biK\ sin(7* + |). 
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Therefore 


Similarly we get 


^2^2 ± dxU 2 = 62 


Thus we conclude 
bi 


= h 


^ sin(7* + 

3 > 



^ sin(7* — 

i) 



sin(7* — 

tt) 

k *2 



at X = . 


at X = ±IX 


at X = ilg . 


= bs 


sin(7* + |) sin(7* — |) sin(7* — tt) ' 
Furthermore, ui{±ll) + U2(±^2) + U3(±^|) = 0 reads as 
Toi sin(7* + ^) + ^1 cos(7* + 

TT TT 

=F 02 sin(7* - -) + 62 cos(7* - -) 

=F 03 sin(7* — tt) + 63 cos(7* — tt) = 0 . 
Altogether, we have to find solutions to the following system 

TT TT 

oi sin(7* + -) + 02 sin(7* - -) + 03 sin(7* - tt) = 0 , 
bi cos(7* + ^) + 62 cos(7* - ^) + 63 cos(7* - tt) = 0 , 


bi- 


= 62- 


■sin(7* + |) sin(7* - |) 
Due to the identities ( 12 . 51 ) and ( 12 . 61 ) we get 


= bs 


sin(7* — tt) 


Oi sin(7* + ^) + 02 sin(7* - ^) + 03 sin(7* - tt) = 0 , 
5 i = 62 = &3 . 

Therefore, in view of the formula ( 12 .6L we obtain 
(01,02,03) e span|(l, 1,1), (0, 


Z(7»-fr 1)}^ (6 i, 62,&3) e span{(l,l,l)}. 


This shows the following lemma: 

Lemma 5.17. Assume K 2 7^ 0 . Then the space of Jacobi functions is a three 
dimensional vector space whose basis consists of 



/cos(k)‘x)\ 

/sin(«;*x)\ 

( ° \ 


C0s(k2x) , = 

sin(fi:2a;) , = 

sin(K^x) 


Vcos(k|x)/ 

Vsin(fi:5x)/ 

\ siiY { Klx ) j 
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We now consider the case K2 = 0 . The general solution of the linearized 
problem is then 

Ui = Oi sin(/t*x) + bi cos(K^a;), U2 = a2X + &2 ) 

uz = as sin(K3x) + 63 cos(k3x) ( = —03 sin(«;*x) + 63 cos(K3a;)) , 

where we used the fact that = —k* in case j* = j. Let us also remind that 
for 7* = ^ we have 


and so <72 ^2 


92 = 


sin(f) 
— 1 . Therefore, 


and I2 = 


sin(f) 

K* 


92^2 ± d^U2 = Ta2 + 


sin(f) 


&2 ± 02 


sin(|) 


^2 


at X = ±^2 ■ 


Taking into account the calculation done previously for ui and us, the condition 
giUi ± dxUi = (72O2 ± dxU2 = 93O3 ± dxUs reads as 


= b 


= ^3 


sin(|) sin(|) sin(-^) 


* / 


— 


-K-i 


-sin(f) 



Therefore, we conclude 61 = 62 = &3- Furthermore, mi(±I*) + ■U2(±/2) + 
■U3(±/3) = 0 reads as 


=Fai sin(|) + bi cos(^) =F 02 


sin(f) 


+ 62 ± 03 sin(|) + 63 cos(^) = 0 . 


Moreover, using the facts that 61 = 62 = bs and cos(^) = — 5, we see that 

bi cos(^) + &2 + cos(^) = 0 . 


In summary, we have to find solutions to the following system 

02 


Ol 


— 03 = 0 , 


&i = 62 = ^3- 

Therefore, 

(oi, 02, 03) G span{(l,0,l),(0,Ki,l)}, {biMM) G span{(l, 1,1)} . 

Lemma 5 . 18 . Assume K2 = 0 . Then the space of Jacobi functions is 3 - 
dimensional and its basis is given by 


= 


/ COs(k3x)\ 


/ sin(K*a;) \ 


( " \ 

1 


0 

, u ( 3 ) = 


\cos{kIx)j 


V— sin(K*x)y 
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5.3.2 The null space N{Aq) is at most five-dimensional 

Next we try to get an upper bound on the dimension of the null space. 

Lemma 5.19. The null space N{Aq) of the linearized operator Aq is at most 
five-dimensional. 

Proof. We have already shown that the space of Jacobi functions is three- 
dimensional. Therefore it is enough to show that there exist at most two inde¬ 
pendent vectors in the null space N{Ao) for which c ^ 0. 

Take any three vector functions € N(Ao) for which the vector 

constants ^ 0, i = 1 , 2 , 3 . Then as 

e { c = (ci, C2, C3) G : Cl -I- C2 -I- C3 = 0} 

which is a two dimensional subspace of there exist scalars a 1,02,03, not all 
zero, such that 

3 3 3 

0 = ^ ^ = T(^ . 

j — 1 i i —1 

Here T is the linear operator defined by the left hand side of ( 15 . 8 |) . Thus we get 
Si=i G J"(r*), in other words, 

3 3 

i=l j=l 

where is a basis of J(r*). This means that the vectors 

are linearly dependent, and this completes the proof. □ 

Indeed, we will prove in Corollarv l 5. 261 below that the dimension of the null 
space is exactly five. 

5.4 Manifold of equilibria 

Our goal in this section is to prove that near p = 0 , which corresponds to F*, 
the set £ of equilibria of the nonlinear system ( 13 . 121 ) creates a smooth manifold 
of dimension 5 . 

5.4.1 Equilibria of the nonlinear system 

Let us first identify the set of equilibria £ of the nonlinear system ( 13 . 121 ) . Ac¬ 
cording to pe£ if and only if for i = 1, 2 ,3 and j = 1, 2,..., 6, 

P G -Bxi( 0 ,i?), 

0 = S'i(p*,p|E*) 

O = 0,(p) 
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on E*. 










Similarly as done in Section 13.21 we can write the first three equations as a 
vector identity on S* and thereby obtain = 0. Thus 

( p & Bxi (0, R), 

p € 6 ^ I 0 = diiPi, pI-s-) on r* , i = 1,2,3, 

[o = 0 ,(/ 9 ) onE*, j = l,..., 6 . 

Taking into account (|E.4I) . the definition of the balance of flux conditions 
05 , 06 and the condition on curvature 04 , by applying the Gauss theorem, we 
see 

r p G Bxj (0, R ), 

p £ £ ^ < Ki(pi, {J po pr)i) are constant, on T* , 

i0,(p)=O onE*, j = 1,2,3,4. 

Therefore, using Lemma |3.1l we conclude: 

£ = I p G Bxi (0, R) : p parameterizes a standard planar double bubble |. 

5.4.2 Level set representation of standard donble bubbles 

Next we represent standard planar double bubbles as a subset of the zero level 
sets of some smooth functions. Let S^iOi), i = 1,2,3, be the corresponding 
circles to standard planar double bubble T = {ri,r 2 ,r 3 }. In other words, 
Ti C Sri{Oi)^ where and Oi are the radius and the center of T^ respectively. 

Lemma 5.20. LetT = DBr^.y^o{0,0). Then 

{ a G : G4a, r, 7 ) = 0 } = Sr,{0,) D T, (i = 1, 2,3), 

where Gi : x (0, 00 ) x (0, —>■ K. are smooth functions defined by 

rsin (7 + |)Gi(cr,r, 7 ) = i sin (7 + |) (^|crp - , 

rsin (7 + |)G 2 (cr,r, 7 ) = i(^sin( 7 - f )|crp - 2 r sin(|)(cr, ( 1 , 0 )) - sin (7 - tt)) 
rsin (7 + f)G 3 (cr,r, 7 ) = i (^sin (7 - 7r)|cr|^ + 2rsin(|)(cr, (1,0)) - sin (7 - f)^ 
with the property that 

Gi+G 2 + G 3 = 0. (5.11) 

The proof is given in Appendix |B] Next let us look at the gradient of G^. 
Lemma 5.21. Let T = 115^,7,0(0,0). Then 

VcrGi{a,r,"f) = ni{a) foraGVi (* = 1,2,3). 
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Proof. It is easy to see that 

a - Oi = - ni{a) 

Ki 


foraer^ (i = l,2,3). 


Using this, we calculate 

rsin(7 + aG2((J,r,^) = sin(7 - |)cr - r sin(|)(l, 0) 

= sin( 7 - f)(CT- 02 ) 

= forCTer2. 

Similarly we get 

rsin(7 + f )V<^Gi(cr,r,7) = for a e Ti. 

rsin(7+ ^^3(0-) for cr e Ta . 

Since rsin(7+ by the identity (| 2 . 5 p we complete the proof. □ 

Furthermore, the following result holds. 

Proposition 5.22. Let T = DBr,-y,o{0,0)- Then 

' drGi{a,r,^) =-1 foraGTi, 

, 9rG2(cr,r,7) =+sin( 7 - 7 r)) fora€T 2 , 

^ drG3{a,r,j) = +5-^^ - sin(7 - |)) for a G Fg. 

Proof. According to Lemma [ 5.201 we have 


Gi{(J,r,'y) = 0 forcrGFi. 

Therefore, differentiating with respect to r in the definitions of functions Gi, we 
observe 

-rsin(7+ |) 9 ^Gi(cr, r, 7) = sin(7 + |)r for a G Fi, 

-r sin(7 + |)9rG2(cr, r, 7) = sin(|)(cr, ( 1 , 0 )) + sin(7 - 7r)r for cr G r2 , 

rsin(7 + |) 9 ^G3(cr, r, 7) = sin(|)(cr, ( 1 , 0 )) - sin(7 - |)r for cr G Fa , 

which finishes the proof. □ 

Similarly we get 

Proposition 5.23. Let T = DBr,-y,o{0,0)- Then 
( 5^Gi(cr,r,7) = 0 

J 5^G2(cr,r,7) = 27ihi(^+f) ( cos (7 - f )|cr|2 - cos(7 - tt)) 

[ 5 .yGa(o-,r, 7 ) = ~ ~ cos(7- f)) 
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for cr G Fi, 
for a G r2, 
for a G Fa. 
















5.4.3 Five-dimensional smooth manifold 

Throughout this section, without loss of generality, we may assume that the 
center of TJ is at the origin of and that the angle 0 * = 0, that is 

Clearly, £ 7^ 0 as p = 0 parameterizes T* = Z?_Br*,7*,o(Oj 0 ). First we 
demonstrate, by applying the implicit function theorem, that every standard 
planar double bubble -Di?r,7,e(ai, 02 ) sufficiently close to T* = Z3i3r*,7*,o(0, 0) 
can be parameterized by some unique vector function p — (pi, P2, Pa) depending 
smoothly on the parameters ai,a2,r,7 and 9 . We continue then to verify that 
the set £ of equilibria is actually a smooth manifold of dimension five. 

Theorem 5.24. Any standard planar double bubble Zli?r,7,e(ai, 02) sufficiently 
close to r*, i.e., (oi, 02, r, 7, 0) G 5 ^( 0 , 0, r*, 7*, 0) for sufficiently small t, can 
be parameterized by some unique smooth vector function p = p(ai, 02, r, 7, 0 ) G 

i?Xi( 0 ,i?). 

Proof. We use the implicit function theorem of Hildebrandt and Graves, see 
Zeidler PPI Theorem 4 .B], with (xo,j/o) = (( 0 , 0 , r*, 7*, 0 ), O), 

X = X Bs^ (r*) X x K , Z = Y, 

y = {p G c'"+“(r*) X c4+“(r*) x c4+“(r;): pi + p 2 + P 3 = o on s*}. 


F:X xY ^ Z, 

((ai,a2,r,7,6»),p) 1-^ {Fi,F 2 ,F 3 ) 


with 

Fi(ai,a2,r,7,6>,p) := pi, pi o {i = 1 , 2 , 3 ). 

Here Gi are the functions stated in Lemma 15.201 and 


'^>^{■,p^,p^opT^){a) = cr + p,(cr)<(cr) +Pi(pr,(cr))r*(cr) for cr G T* , 
where p = J 7 p on E*. Furthermore, 


Qe = 


cos 9 sin 9 
— sin 9 cos 9 


TffV = V — a 


are the clockwise rotation matrix and the translation operator respectively. 
Indeed, the image of the function F lies in Z = Y, that is 


Fi T F2 T F3 — 0 on S*. 


( 5 . 12 ) 
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To see this, note that for cr S E*, 

Opri)(CT) = ^' 2 (-,P 2 ,Ai 2 Opr 2 )(cr) = 4'3(-, P3, ^*3 O Pr 3 )(CT) , 

by Lemma IXTl This together with the identity (15.111) proves (15.121) . 

Moreover, since according to Lemma [5.201 we have 

Fi{xo,yo){<7) = Fi((0,0, r*,7*,0),0)(cr) = Gi{a,r*,j*) = 0 for cr G T*. 

Thus F(xo,yo) = 0. Now let us compute the derivative dpF(xo,yo)' 

dpFi(xo,yo)(v)(cr) = V„G,(cr, r*, 7 *) • (uj<(cr) + (Jv(pi^(a))).T*(a)) = vi , 

where we used Lemma 15.211 Thus 

dpFixo,yo) = I ■ (5.13) 

Furthermore, F is a smooth map on a neighborhood of {xo,yo)- 

Therefore, according to the implicit function theorem, there exist neighbor¬ 
hoods U = B^{xo) of xq and V = Bxi (0, R) of yo = 0 and a smooth function 

p:U —>V 

(ai, 02 , r, 7 , 0) H> p(ai, 02 , r, 7 , 0 ) , 

such that p(xo) = 0 and for every (oi, 02 , r, 7 , 0 ) G i?g(0, 0, r*, 0) we have 

F((ai, 02 , r, 7 ,0), p(oi, 02 , r, 7 , 0 )) = 0 ■ (5.14) 

Moreover if {x,y) G U x V and F{x, y) = 0 then y = p(x). 

We now claim that Tp = {rp^,rp^,rp^} parameterized by the function p = 
p(ai, 02 , r, 7 , 0) is the standard planar double bubble DBr,-yfi{ai,a 2 )- To see 
this, note 

Fj((oi,02,r,7,6»),p(oi,02,r,7,6»)) = 0 

Gi(^(5eTs4'i(-,pi,/riOprJ,r,7^ = 0 
QelffTp. C S'r-. (Oi) by Lemma 15.201 

Therefore, since Lemma ld.ll guaranties that the curves Tp^, Tp^, Tp^ meet at their 
boundaries, we end up with two choices: Either Tp^ = Tj, where T = {Ti, r 2 , Ta} 
is a standard double bubble DBr^-y^eicii-, 0,2) or Tp^ is the complementary part 
of Ti in SniOi). But the latter can not happen since the norm of p is small. 
Hence 

^^(ai,02,r,7,0) F)Bj-^^ g{^Cl\^ O2) , 

as required. □ 

Theorem 5 . 25 . The set of equilibria £ is in a neighborhood of zero a C^- 
manifold in Xi of dimension 5. 
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Proof. Remind that we have shown 

f n [/ = |p G Bxi (0, R) : p parameterizes a standard planar double bubblej fl U 
= |p(ai,a2,r,7,6») : (m, 02, r, 7,6») & U = ^^(0,0, r*, 7*, 0) |, 
where the function 


p: u ^Xi = C'‘+“(rt) X C'^+“(r*) x c^+“(r5) 

(ai,a2,r,7,6>) i-)> p (ai,a2,r,7,6») 

is smooth, in particular and p{U) = f, p{xq) = p(0,0, r*, 7 *, 0) = 0. 

Therefore, it is left to check that the rank of p'{xo) equals five (see the 
definition of a manifold on page usi. To do this, we differentiate (15.1411 with 
respect to i G {oi, 02 , r, 7 , 0 } and evaluate at xq to get 


d,F{xo, 0) + dpF{xo,0)di,p{xo) = 0. 


Therefore, (15.131) gives 


d,p{xo) = -d,F{xo,0) (t G {ai,a2,r,-f,9}). 


We now calculate 


da^F,{xo,0) = VaGi{a,r* ■ (-1,0) = n*{a) ■ (-1,0) = cos{k*x) , 

where we used the fact n*{a) = —(cos(K*a;),sin(K*a;)), i = 1,2,3. Thus 

daiP{xo) = (cos(K*a;),cos(«;2a;),cos(K3a;)). 

Similarly, we get da^p{xQ) = (sin(K 3 a;),sin(K 2 a;),sin(K 3 a;)). 

Next we calculate 

(xo, 0) = (a, r*, 7*) • ( [ 1 ] • a) 

= <(c^) ■ (^^ = n*[a) ■ { - ^n*(cr) + O*)^ 

= <(a)-0*^ 


and so 


9oPixo) = 


sin(K^a:) 
\ sin(K|a:) 


We now compute the derivative drF{xo,0) = drG{a,r*,'y*). According to 
Proposition 15.221 


5rG2(cr,r*,7*) = 


1 

sin( 7 * + |) 



tTi + sin( 7 * 
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First we consider the case ^ 0. Employing the arc-length parameterization 
of r 2 derived in Proposition 1C . 1 1 we obtain 

E, ^ - /sin(f) 

drG 2 [a,r ,7 ) =-'-^ 


V 


sin( 7 * -h f) 
kI sin(^) 


K* sin( 7 * -t f) 

sin(f) 


cos(K 2 a;) — 


(Ti -I- sin( 7 * — tt) 

1 / sin^(|) 


( 


sin( 7 * -I- |) Vsin( 7 * - |) 


-|- sin( 7 * — TT 


, , sin(7*-hf) 

COSiKnX) -—-^ 

^ 2 ' sin( 7 * - |) 


sin( 7 * - |) 

where we applied the formula sin^(a:) — sin^(j/) = sin(a; -|- y) sin(a; — y). 

A similar argument works for drG^{a,r *Altogether we derive in case 

TVf, 

/ 1 \ 


drp{xo) = 


sin(f) ^ ^ sin(7* + §) 


sin(7* - I) 

sin(f) 


— cos(k 2 x) + 


\ sin( 7 * — tt) 


sin(7* - I) 

sin( 7 *-hf) 
cos(kSx) + . , ^ ^ 

^ ’ sm( 7 * - tt) / 


Next we consider the case K 2 = 0- We calculate 


5,G2(^T,rM) = - 


1 /sin(|)r* 

sin(f) V r* y 




sin 


3 V 


Therefore, we derive in case k% = 0, i.e., when xq = (0,0, r*, 0), 

( 1 

drPixo) — I 2 

y— cos(kJx) — \j 

Finally let us calculate i9^F(xo,0). We have djF(xo,0) = 9^G(cr,r*, 7 *). 
We first consider the case K 2 0: Employing the arc length parameterization 
of F* to the formulas derived in Proposition 15 . 231 we derive in case K 2 0 that 

/ 0 

9^p(xo) = 02 COs(k2x) -I- &2 

\a3 COS^K^x) + 

for some constants ai, bi (see the Appendix for the explicit form of the constants). 
This immediately implies that d^p{xo) is independent from the other elements 
of p'{Xo). 

However, we give the explicit formula in case K 2 = 0- Using Proposition [523] 
we see 


a^G 2 (a,rM) = 


1 


2r* sin(^) 




/I 


+ 77- 


3 1 


sin(f)V2 8(«)E) 
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d^Gz{(T,r 


* TT 'i 

3^ 


(cos(-f)|ap-(r*f) 


2r* sin(^) 

2 r*sin(^) 2 

-+ \ cos{kIx)) . 


Ki sin(-|^) 


In summary, we have proved that the rank of p'{xq) is equal to five and we 
have shown that the set of equilibria f is a C'^-manifold in Xi of dimension five. 
Moreover 

where 


z;(i) = 


/cos(Kja;)' 

cos(k 22 ;) 

\cos{kIx)^ 


z;( 2 ) = 


/sin(Kja;)' 

sin(K 2 a;) I , = 

\sm{K^x)^ 


I 


0 


sin(7*) 
sin( 


(" .-D sin(K^a:) 
\ sin(/C 3 a:) 


^;(4) = 


( ' \ 

— sin(-^) / * \ , sin{ 7 * + ^) 

. , cos(/cSx) + . ; . 

“(7 - 3 ) ^ ^ ’ sin( 7 *- 3 ) 

^™(f) pns('K-*Ti + + / 

\sin(7*-7r) COb[K^X) + sin(7*-7r) / 




/ 0 
02 C0s(k2x) + &2 I . 
\a3 cos(K5a:) +63^ 


Although are continuous in particular at K 2 = 0, for convenience we state 


them in case = 0: 


^(i) = 


^C 0 s(k 3 x)' 
1 

, cos(K 3 a;), 

/ 




= 


1 


y— COs(K 3 a:) — 1 ^ 


sin(K*x) 

0 

sin(K*x) 


= 



/ 


\ .lsin{^) (l^OsiKlx) + l)J 


□ 


5.4.4 Geometric interpretation of the null space 

As an immediate corollary of Theorem 15.251 we get 

Corollary 5.26. The null space N{Aq) is five dimensional. Furthermore, 

To£ = N{Ao) . 
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Proof. It always holds 


To£ C N{Ao), 

see [TJ equation (2.8)]. Thus, according to Theorem 15.251 and Lemma [5.191 
5 = dim(Tof) < dim(A^(Ao)) < 5 . 


It follows that dim(iV(ylo)) = 5 and moreover Tq£ = N{Ao). □ 


Variations preserving areas and curvatures 

We easily see, using formula (12.5|) . that 



In other words, 


j{r*) c j-(r*). 


(5.15) 


By Lemma 15.71 each of the (i = 1,2,3) corresponds to a first variation of 
r* which preserves the areas, and the curvatures to first order. Indeed, we have 
demonstrated in the proof of Theorem l5.25l that correspond to the 

first variations of the double bubble T* associated with translation along x-axis, 
translation along j/-axis and rotation around the center of T^, respectively. 


Variations not preserving areas and curvatures 

It is shown in the proof of Theorem 15.251 that corresponds to the first 
variations of the double bubble T* associated with uniform scaling (with the 
scale factor ^). Let Ai{r) denote the area of the regions Ri{r) corresponding 
to the double bubble ZlBr, 7 *,e*(a*, aj). Then (see equation (3.1) in [12] ') 

drAi= f — f drA2= f 2 — j (5.16) 

according to Lemma ID. 11 (ii). 

Again remember from the proof of Theorem 15.251 that corresponds to 
the first variation of T* with respect to the angle 7 , that is w.r.t. the curvature 
ratio. Similarly we denote by Ai{"f) the area of the regions Ri{'j) corresponding 
to the double bubble DBr*, 7 .e*(oi,oj). Then 

d^Ai = f - f >0, d^A2 = f v^^^2 - f < 0 (5.17) 

JT* *^^2 '^^2 *^^3 
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according to Lemma ID. 21 (ii). 
We now define the matrix 


/drAi 

d.yAi\ 

= 

(ir-j 

f V^'^^2 

f — f 

Jri Jr- 


d~fA2j 


■ j 

f 



Lemma 5.27. The matrix D is invertible for each 0 < 7 * < 

Proof. Let us calculate its determinant. Inequalities (15.161) and (15.171) imply 
det D = drAidjA 2 — dyAidrA 2 < 0 

Now as the determinant of the matrix D is strictly negative, we conclude that 
the matrix D is for each 0 < 7 * < ^ invertible, □ 

As a further result of Lemma [D! T] and [D] 2 ] (ii), we get ^ J^(r*). 

Therefore, we conclude from Lemma ISJl that the corresponding variations do 
not preserve areas to first order. Indeed we will show below in Lemma lA. 21 that 

l{u,u)<0 ioT u = 

In addition they do not preserve the curvatures to first order too as the constant 
vectors c{vA')) and are nonzero. 

5.5 Semi-simplicity 

We need to show two small propositions. The first one is stated and proved in 
the proof of Lemma 3.8 in |12| . 

Proposition 5.28. If u G J^(r*) satisfies I{u,u) = 0, then 

/(u,u )=0 VueJ^(r*). 

Proof. According to Corollarv l5.12l I{v + tu,v + tu) > 0 for all v G J^(r*) and 
t gR. Therefore 

I{v + tu,v + tu) = I (u, v) + 2tl{u, v) + t^I{u, u) 

= /(u, v) + 2tl{u, v). 

This forces I{u, v) = 0 as t can take arbitrary negative values. □ 

Proposition 5.29. Let z G R{Ao). Then there exists u G P(r*) (h L)(Ao) such 
that Au = z. 
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Proof. Clearly, there exists u G D{Aq) such that Au = z. The actual task is to 
find two constants a{u), f3{u) such that 

u :=u + + /3{u)v^^^ 


satisfies 


-'ll -'i 2 -'is 

(This will finish the proof since € N(Ao) implies Au = Au = z.) To 

do so, let us recast this integral constraint into the matrix form 



where the matrix D is given above. The invertibility of this matrix proved in 
Lemma [5 . 271 finishes the proof. □ 


We are now ready to prove: 

Lemma 5.30. The eigenvalue 0 of Aq is semi-simple. 

Proof. Since the operator Aq has a compact resolvent, the semi-simplicity condi¬ 
tion is equivalent to the condition that N{Aq) = N{Aq) (use the spectral theory 
of compact operators, e.g. see [4, Section 9.9]). In other words, it suffices to 
check that 

i?(Ao)nA^(Ao) = {0}. 

To prove this, let 2; € R{Aq) fl N{Aq) ( C D(Ao)). According to Proposition 
15.291 there exists u G T>(Ao) fl J^(r*) such that Au = z. From this, exactly as 
done in Section [5.21 we derive the identity 



Vr*(Ar*Ui -f (K*)^Ui)j^ds = I(z,u), 


(5.18) 


where we used only the facts that u,zG D(Aq). 

Moreover, similarly as before, an integration and application of the diver¬ 
gence theorem using the fact that u G D{Aq) gives 



and so 0 e J^(r*). 

Now since z € iV(Ao) fl J^(r*), Corollarv l5. 161 tells us I{z,z) = 0. Therefore, 
according to Proposition 15.281 


I{z, u) = 0 

as u G J^(r*). In view of the identity (I5.18|) . we obtain u G N{Aq). Conse¬ 
quently z = Au = 0, which finishes the proof. □ 
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Remark 5.31. The main ingredient in the proof is the positivity of the bilinear 
form, i.e., the variational stability of the stationary solution. 

6 Standard planar double bubbles are stable 

Summing up, we have shown that all the hypotheses of Theorem ll. H are satisfied. 
Thereby applying Theorem 14.11 we conclude: 

Theorem 6.1. Let T* be a standard planar double bubble. Then p* = 0 is a sta¬ 
ble equilibrium of (13.121) in Xi = R^). Moreover, if is sufficiently 

close to p* = 0 in Xi and satisfies the corresponding compatibility conditions 
(031), then the problem (13.121) has a unique solution 

peC^+^’^"^+°‘{[0,oo) xH, 

and approaches some p°° G £, parameterizing some standard planar double 
bubble, exponentially fast in Xi as t ^ oo. 

In this sense, the standard planar double bubble T* is stable under the 
surface diffusion flow. In addition, every planar double bubble that starts suf¬ 
ficiently close to r* and satisfies the angle, curvature, balance of flux condition 
and the condition on the Laplace of the curvatures, see (12.2L at t = 0 exists 
globally and converges to some standard planar double bubble, enclosing the 
same areas as its initial data, at an exponential rate as t —>■ c». We illustrate 
this result in Figure [TJ 

6.1 General area preserving gradient flows 

It is to be expected that for any sufficiently smooth area preserving gradient 
flow 

V = —grad^(p) Length 
one obtains the following identity 

lbllw(r*) = , (6.1) 

where z := (5(grad.j^(p)Length) (u). Here 'H(r) denotes a (pre-)Hilbert manifold 
with some area constraints. 

In particular, if u is a eigenvector of the operator (5(grad.j^(p)Length) with 
respect to the eigenvalue A, then we get 

= XI(u, u). (6.2) 

w(r*) 

Comparing the identifies (16.21) and (j6.1l) with the identities (j5.5l) and (j5.18|) 
respectively, we expect that our approach can be used for other area preserving 
gradient flows. Therefore we conjecture that 

Conjecture 6.2. Standard planar double bubbles are stable under sufficiently 
smooth area preserving gradient flows. 

It would be desirable to analyze the problem systematically. 


(5 (grad^(P) Length) (u) 
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Appendix A More about the bilinear form /(, ) 

Lemma A.l. Within the class of functions u satisfying the linearized angle 
condition, we have 

{ u : i{u, w) = 0} n Tir*) = N{Ao) n T{r*). 

Proof. We first assume u G J^(r*) such that I{u,u) = 0. Then by Lemma 3.8 
in m and the fact that u satisfies the linearized angle condition, we conclude 
that u G N(Ao). The converse statement is Corollary 15.161 □ 

Note that we have already shown in Section 15.4.41 that 

N(Ao)nP(r*) =J(r*) =span{v(^\v(^\v(^^}. 

On the other hand we obtain: 

Lemma A. 2. For the bilinear form I it holds 

l{u,u)<0 for u = (0 < 7 * < ^). 

Proof According to Lemma 15.151 

3 

I{u, u) = -y^ Cjju) 

Now assertion (iii) in Lemma [DN] and Lemma [0.21 proves the lemma. □ 



Appendix B The proof of Lemma 15.20 


Consider the standard planar double bubble T = Zli?f., 7 ,o( 0 ,0). That is the 
left circular arc Ti has radius ri = r centered at Oi = ( 0 , 0 ) and all the other 
centers also lie on the x-axis, for some r>0, 0 < 7 < see Figure IH] 

It follows directly from the law of sines that in case 7 ^ f 


O 2 

O 3 


Vsin( 7 -f) ’ V’ 

C^r,Q] 

V sin( 7 ) ’ / 


\ 

sin(^ - 7 ) 


sin (7 - f) 


?'3 = 


sin(^ - 7 ) 
- - - 1 

sin( 7 ) 


Therefore, for a = [<Ji,<J 2 ) G r 2 , 7 we have 


0 = 


(s 


sin(^) 
sin (7 - f) 


sin(^ - 7) \2 


r 0^ - ^ -—r) 

' ) Vsin(7-?) / 


= lal^-2a 
2 

sin (7 - f) 


sin(^) 


sin (7 - f) 

sin^(¥)-sin"(f-7) 


( " l. r,0) + ^ 1.^ - -r^ 

\s 


.sin (7 - f) 


rsin( 7 + §)G2((7, r,j), 


sin^( 7 -f) 


39 





























Figure 6 : Standard planar double bubble F = I?_Br, 7 ,o( 0 ,0) 

where we applied the formula sin^ x — sin^ y = sin(a; + y) sin(a: — y). 
Similarly, for tr S Fa we obtain 


0 = 


V sin(7) ’ / V sin(7) 


sin(f-7) 


sin(7 — tt) 


rsin( 7 + f)G' 3 (cr, r, 7 ) 


and obviously for a G Ti we have 


0 = lo-r-r^ = 


sin (7 + f) 


rsin (7 + f )Gi(cr, r,j). 


Furthermore, we see for cr G F 2 , 7 = | that 0 = | — cti = G 2 (cr, r, |). 

Finally, the identity (12.61) easily verifies (15.111) . This finishes the proof as the 
coefficients appearing above are all nonzero and well-defined. 


Appendix C Arc-length parameterization of T* 

Proposition C.l. An arc-length parameterization ofT*, z = 1,2,3 is given as 
follows: For 7* ^ 


1 


(cos(K 3 a;), sin(/cja:)) 


for CT G F^ 


(o-i,o- 2 ) = cr = < 


(s 


.y,* _|_ J_ (;.os(/c*a:) , ^ sin(«;2a:) ) foraGV^, 


sin( 7 * - f) K* 


— -r* + ^ cos(K 3 a;) , ^ sin(K 3 a;)) /or cr G F 3 , 


Asin(7*-7r) k* 

Moreover this arc-length parameterization is continuous at "f* = ^ and in par¬ 
ticular cr = (^ , x) for cr G F 2 , 7 * = ^. 
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Proof. We give the proof for . We observe 


(cTi, 0 - 2 ) = cr = O 2 - = O: 


■ (cos(K2a;), sin(«:2a;)) 


^2 

sin(f) 


sin( 7 * - |) 


^ ^ cos { k 2 x ) , ^ sin(K^a;) 


for (T G Fn 


The proof of the continuity can be done using the identity (12.51) and the L’Hopital’s 
rule. □ 


Appendix D The signs of the integrals 


Lemma D.l. Let 0 < 7 * < Then 


(i) f f u^"‘^ 2 < 0 , f 

t/r? t/To t/ro 


(ii) / U^l-/ U^"^ 2 > 0 , / V^^>2- I V^^>3>0, 

3 


IT. 

,( 4 ) 




,(4). 


(iii) ^Ci{v^) f > 0 . 

i=i 


Proof. In order to easily see the general strategy of the proof, let us first verify 
the assertions for 7 * = ^: 


II Al 

= l = 2ll>0, 

J-l* 


= (k*)^ > 0 , 
02 ( 1 ^^"^^) = 0 , 

[ = — f 1 + cos(Kja;) < 0 , C3(u*'‘*^) = — (wg)^ < 0 . 

J-lt J-lt 


f-l 


'-1% 




Therefore, 


>0 ( 7 * = f)- 


Next we calculate 

f ^ =-l* + 21* + f cos{kIx) = -1*2 + 211 + 

Jt‘ J-q 

sin(f) ^ sin(f) 


^,( 4)2 - 


+ 2/t-2- 


= 2^1 +r* sin(f) > 0 , 
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where we have used the facts that 


/* 7* /* 

n — ^3’ ^2 — 


sin(f) 


^*J* _ 27r 

, — “3" 


(7* = f). 


Obviously Jp. v^'^\ — /p. > 0 which completes the proof of assertions (i)- 

(iii) in case 7 * = |. 

Assume now 7 * f • Then we calculate /p. = /p. 1 = 21* > 0 and 


i r uw.= 1 

2j_i, 2j_^Asi 


- 1 % Vsin(7* - f) 


f) ^ ^ , sm(7* + f)' 

^ cos(K2a;) + 3-7—-^ 


sm(7* - f). 


-sin(f) sin(K^l^) sin( 7 * + f) ^« ^ sin(f) sin( 7 * + f) ^ 

sm(7* - I) sin(7* - f) ^ ^2 sin(7* - f) 

, /-sin(|) ^ sin( 7 * + f)x _ 


. 3 J sm^7 + 3 J X ^ 

f) ^ sin(7*-f)J • 2 /( 7 ), 


where we used the fact that K 2 I 2 = —{l* ~ f)- Similarly 




sin(f) sin( 7 * + f)\ ^ 

, , ^ \ + —tV^) =: 21 * 3 7* . 
( 7 * — tt) sm( 7 * — tt) / 


Obviously the function 3 is negative on (0, ^). Taking into account the fact 
that sinlx) < a; for 0 < a; < tt, it is easy to check that f( 0 ) < 0 , f' < 0 and so 
/ < 0 on (0, too. Thus 

= 2l*2f{r) < 0 , = 213*3(7*) < 0. 

J-q J-q 

Assertion (i) follows. 

Similar argument shows that 3 ' > 0 and so /' — 3 ' <0 on (0, ^). This 
together with (/ — g){^) = 0 implies / — 3 > 0. Observe further that 


sin( 7 * — tt) ^ sin(—sin(|) sin( 7 * — |) 
( 7 *-^) (7*-f) 


on ( 0 ,f) 


as the function 2 H^£) ig strictly increasing and decreasing on intervals {—tt, 0 ) 
and (0, tt) respectively. Thus we conclude 


Is ^ ( 7 * - tt) ^ (7*-7r) sin(7*-f) 

^2 (7*-f)«3 sin(7*-7r) ( 7 * - f) 


(0<7*<f). 


We are now ready to estimate 




f " = 212/(7*) - 213*3(7*) 

J-li 


>212(/(7*) - 5(7*)) > 0 . 
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Moreover, Jp, — /p, > 0 by assertion (i). This proves assertion (ii). 

To prove assertion (iii) we observe: 

^ (4)^ , *^2sin(7* + f) , * , . *^2sin(7* + f) * * 

C 2 {v^ 0 = (^2) ^ = K1K2 , C 3 (r;^ >) = (K3) ^ = K1K3 , 

sin( 7 * — sin( 7 * — n) 

and = (k*)^. Therefore, taking into account that k* < 0 

3 „ 

^ aiv^) J ^ = 211(^1)^ + 2i;Kl4f{^*) + 2i;Kl4g{j*) 

> 2Z2*«t«*/(7*) + 2;3*<«*/(7*) = 2Klf{j*){l*24 + 1*34) 
= 2 Kt/( 7 *)(-( 7 * - I) - ( 7 * - tt)) = 4«t/(7*)(-7* + f) 


> 0 . 


Lemma D.2. Let 0 < 7 * < Then 


(i) f v^^\ = 0, j u^®^ 2 < 0 , f v^^^3>0, 

t/p? 


(ii) 




> 0 , 


/r? 


[ - [ V^^'>3 < 0, 

Jr* Jr* 


^ r 

(iii) ^Ci{v^) / > 0 . 

Jr* 


□ 


Proof. Let us first consider the case 7 * = f • Then 

f ^(5)3=0, ci(vi^i)= 0 . 


i-it 

pf-o 




rl* 


l-l* 

pit 


sin(f) 


' <». c4.'7 = a^<0, 




’-K 


8 {Kty 


= 


J-l* 

Therefore, 


/_) i c<,s(k:.) + 1 > 0, «(."■>) = > 0 . 


SZ' 

2=1 


.(5) 


*>0 (7* = i)- 


r* 


Assertion (ii) is an immediate consequence of assertion (i). This proves (i)-(iii) 
in case 7 * = f • 
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Now assume 7 * ^ Clearly Jp, = 0 as = 0. Next we compute 


kP = 


+ 


cos(K 2 a;)- 


, Op* sin(f) 

(^*)2 ' sin 2 ( 7 *-f)^ ^ sin( 7 *-|) 

sin"(7*+ 2 , sin"(f) .„*^2 !in(i)sin(7* + f) 


sin" ( 7 * - f) 


{ry + 


-{r*Y-2{r*Y 


3 , sin"( 7 *-f) 

■'-'aM- foraer;. 


sin"( 7 *- f) 


cos(«: 2 a;) 


> .^*^2^ sin(7* + f)-sin(f) 2 


sin( 7 * - 


Therefore, 


-/ « 


2=1 djG2{(J,r*,j*) 

= r* sin(V + f) /^cos(7* - f )|a|2 - (r*)2cos(7* - tt) 

j 


/r^ 


where 


f{x) := cos(x - |)( 


^,^^sin(x + f) - sin (|),^2 


) + cos(x). 


sin(x - I) 

It is not hard to show that the function / is strictly decreasing. Together with 
the fact that this function vanishes at 7 * = ^ we conclude / > 0 and so 
/„.<»!. < 0 . A similar proof works for This completes the proof of (i). 

The statement (ii) is an immediate consequence of assertion (i). Similarly 
you can check that C 2 (w^®^) < 0 and 03 ( 1 ;^®^) > 0 which easily gives (iii). □ 


Appendix E Deriving the parabolic system 

For the normal velocity Vi of ri(t) := Tp^^^. (t) we obtain with the convention 

em 

= i, Rda'^i)dtpi{a,t) + a G T*, 

where the unit normal rii of ri(t) := rp.^p;(f) is given by 


ni{a,t) = —Rdc^i{a,t) 

= ^{Rda'^^ + Rd.ui'^idaPi{(J,t)) cr G r* . 
Ji 


(E.l) 


Here 


Ji = Ji{cr,pi,pi) := \dc^i\ = \/\dc,'^i\^ + 2{da'^i,dw'^i)dc,pi + \dw'i'i\'^{daPiY, 
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and R denotes the anti-clockwise rotation by 7r/2. Next computing the curvature 
Ki(= K^{(j,pi,Hi)) of rj(t) := we get 




1 

J" 

1 

J .3 1 


(E.2) 




Therefore, the surface diffusion flow equations can be reformulated as 

dtPz = ai{(T, Pi, pi)l^{a, p,, pi)Ki{a, p,, Pi) + bi{a, Pi, pi)dtp^, (E.3) 


where 


ai{cj. Pi, Pi) 


J,{a, pi,p{) / ^ _ 1 _\ 

(d^u'l'i,Rda'i'z) ^ (n*(a),ni(a,t)) / ’ 


bi{a,pi,pi) := 


{dr'i’i,Rda^i) + {dr'^i,Rd.u,^i)d^Pi / {T*{a),ni{a,t)) 


-{d^'^^,Rd^^i) 


( 


{n*{a),ni{a,t)) 


A{a, pi,pi)v 


zd. 


Ji)^^Pi-)Pi) ^ Ji{,^i Pi-: pi) 



Note that we have omitted the projection pr^ in the functions pi and the term 
{a, pi{a,t), pi{pTi((T),t)) in du'^i with u G {a,w,p'\ to shorten the formulas. 
Furthermore note 


6*Ip.^o = -(t;,±<) = 0, a,|p,^o = l. (E.4) 

We will now make use of the linear dependency (13.91) to derive from the equations 
(IE.31) evolution equations solely for the functions pi. For this, let us rewrite (IE.3|) 
into 

dtPt^^t{.Pi,p\T.-)+ '^i{pi,p\Y.‘)dt{jpopv^)^ inE*, (E.5) 

where for cr G E* 

?*( p *, p | e *)( o ') = ai{a,p^, {Jp\^>^)i)A{a, pi, {Jp\^*)i)Ki{a, pi, iJp\^*)i), 
‘Bz(Pi,p|E*)(o-) = h{a,pi, {Jp\j^*)i), 

and where we used the linear dependency (13.91) . By writing (IE.51) as a vector 
identity on E* we get 

dtP = ^{p-,p\Y.-)+^{p,p\Y.-)J{dtp) on E*, (E.6) 
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where we employed the following notations 


S'(p,p|s*)(o-) := (5 '*(p»,p|e*)(o')) for cr e E* , 

V / 2=1,2,3 

‘B(p,p|e*)(o-) := diag((«8i(p*,p|E.)(o-))i=i,2,3) foi' o’ e E* . 

We rearrange to find 

{l-^{p,p\^,)j)dtp = ^{p,p\^^) onE*. (E.7) 

Consequently we get 

onE*. 

According to (IE.4I) . in some neighborhood of p = 0 in C^(r*) the inverse (/ — 

^{p, p\y.*)J) ^ exists. Inserting the above equation into the equation (IE.61) we 
can finally reformulate the surface diffusion flow equations 

Vi =-Ar^Ki onri(t) 

as a system of the evolution equations for functions pi defined on fixed domains 
r* (or equivalently on 

dtPi = ^z{pi,p\Y.-) +*Bi(pi,p|E*)({^(d- ®(p,p|e*)^) ^S'(p,p|e*)} °P^*) ■■ 

Finally, we rewrite the boundary conditions at cr G E* as 

®i(p)(c^) := Pi (o') + P2(cr) + p3(cr) = 0 , 

®2 (p) (cr) := (ni (cr), n 2 (cr)) - cos ^ 

= + dyj'HidaPl) , ■^idcr^2 + dy,'^2 daP 2 )) - COS ^ = 0 , 

Jl J 2 

®3(p)(cr) := (n2(cr), n3(cr)) - COS 

= (^(5<t^'2 + dy,^2 daP2) , ^(9a4'3 + 9i«^'3 d^Pi)) - COS ^ = 0 , 

3 

04 (p)(cr) :=^Ki(cr,pj,(J^p|E*)0 =0> 

2 = 1 

® 5 (p)(cr) := {J p\t.*)i) - ^da{K2{(J, p2, {J p\y.-)2) = 0 , 

Jl J 2 

® 6 (p)(cr) := -j-da{K2{cr,P2, {Jp\^-)2) - —da{K 3 {a,p 3 , (J^pIeOs) = 0 . 

J 2 J 3 

We emphasize that the operators ©i {i = 1,..., 5) are purely local due to the 
fact that the projections pr^ act as the identity on their image E*. 
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